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State selective protocols, like entanglement purification, lead to an essentially non-linear quantum 
evolution, unusual in naturally occurring quantum processes. Sensitivity to initial states in quan¬ 
tum systems, stemming from such non-linear dynamics, is a promising perspective for applications. 
Here we demonstrate that chaotic behaviour is a rather generic feature in state selective protocols: 
exponential sensitivity can exist for all initial states in an experimentally realisable optical scheme. 
Moreover, any complex rational polynomial map, including the example of the Mandelbrot set, can 
be directly realised. In state selective protocols, one needs an ensemble of initial states, the size of 
which decreases with each iteration. We prove that exponential sensitivity to initial states in any 
quantum system have to be related to downsizing the initial ensemble also exponentially. Our results 
show that magnifying initial differences of quantum states (a Schrodinger microscope) is possible, 
however, there is a strict bound on the number of copies needed. 


Introduction 

Quantum technology progresses at a fast pace. Prepa¬ 
ration, control and measurement of coherent quantum 
systems [1] became possible on an unprecedented level 
leading to a wealth of proposals of applications ranging 
from quantum information processing to high precision 
measurements and sensors. In these protocols, increas¬ 
ingly sophisticated sequences of coherent evolution, mea¬ 
surement and post-selection are applied in order to con¬ 
trol the state of quantum systems. Dynamics achieved 
by state selective protocols was proven essential for a 
large number of quantum information protocols [H 13] 
and quantum communication [3| . Prominent examples of 
probabilistic protocols are the KLM scheme [5] for linear 
optical quantum gates or the entanglement purification 
protocols mm employing measurement and selection in 
order to increase the entanglement between subsystems. 

Manipulation by measurement and selection breaks the 
linearity of quantum mechanics, thereby broadens the 
possibilities for quantum evolution [HHII]. In contrast 
in the well established field of quantum chaos [T3| one 
studies the signatures of chaos in closed quantum sys¬ 
tems with linear evolution. However, the essential non¬ 
linearity of an iterated, state selective protocol can result 
in truly chaotic behaviour, showing exponential sensitiv¬ 
ity to initial conditions [ISHTS]. So far sensitivity has 
been proved only for a tiny fractal subset of initial states 
with zero measure. In this article we demonstrate that 
exponential sensitivity can exist for all initial states in 
an experimentally realisable optical scheme. Moreover, 
we show that any complex rational polynomial map, in¬ 


cluding the example of the Mandelbrot maps m, can be 
directly realised using state selective protocols bringing 
a whole new class of quantum protocols to life. 

From a fundamental point of view, one can search for 
the most general evolution for a quantum system. A very 
general dynamics is sometimes imagined as a system to¬ 
gether with one or more ancillas and allowing for both 
unitary evolution and non-selective measurements on the 
complete arrangement. The evolution reduced for the 
system only is called a quantum channel. When talking 
about quantum states in practice, it is unavoidable to be 
able to repeat experiments on an ensemble of identically 
prepared initial states, in order to uncover the underly¬ 
ing probabilistic laws. This ensemble view of quantum 
states allows for the following trick when designing the 
most general dynamics for a given initial state. Let us, 
for example, consider systems from the ensemble pair¬ 
wise and let them interact with each other. After the 
interaction one can perform a measurement on one of 
the pairs and then discard the measured member of the 
pair. In case of selective measurement, one may also 
discard the unmeasured member of the pair, depending 
on the measurement result. The resulting ensemble will 
be reduced in size, but some of its properties may be 
changed in a benehcial way, e.g. entanglement between 
subsystems. The above procedure goes beyond the usual 
notion of quantum channels, in the following sense. The 
initial step of the procedure, namely taking the systems 
pairwise, can be viewed as splitting the original ensemble 
into two parts, and employing one part as an ancilla. In 
other words, the state of the ancilla will be dependent on 
the state of the system. The state dependent ancilla lies 
at the heart of the non-linearity of the process. 
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Results 

A linear optical experimantal scheme implement¬ 
ing a family of non-linear maps. We propose here 
a simple experimental setup which implements a non¬ 
linear process exhibiting exponential sensitivity to the 
initial state. Our scheme is inspired by an experimen¬ 
tally tractable entanglement purification protocol HZIIIH] 
and uses only linear optical elements. During iterations 
we form pairs of photonic qubits from an ensemble of 
identically prepared photons and apply a post-selective 
transformation on the pairs by measuring the polariza¬ 
tion of one photon and keeping or throwing away the 
other photon depending on the measurement result. The 
post selection induces a non-linear, deterministic trans¬ 
formation on the remaining photons, therefore the kept 
photons remain identically prepared. 



(a) PBS {h)\H\\H), {c)\V\\V), id)\H\\V}, {e)\V\\H), 





(g) Simplified scheme 


(h) Two iterations 


Figure l:The proposed experimental setup, (a) A Po¬ 
larizing beam splitter (PBS) with two spatial input/output 
modes. |(b)|(e)| The effect of the PBS acting on the four pos¬ 
sible two photon input states regarding polarizations. |(f)| 
A post-selective linear optical scheme inducing a non-linear 
transformation (g) and its simplified version utilising a polar¬ 
izer reducing the success probability by 1/2. (h) A two level 
scheme amended by a unitary transformation U acting on the 
polarization state of the photons. We consider a run of this 
experimental setup successful if all the detectors click. This 
condition introduces the post-selection to the system. 


Let us denote the horizontal and vertical polarization 
states \H) and \V) for our photonic qubits. The key ele¬ 
ment of our scheme is the polarizing beam splitter (PBS). 
When two photons arrive at the same time but from dif¬ 
ferent spatial input modes this linear optical element in¬ 


troduces entanglement between the spatial modes and 
the polarization degrees of freedom, see Fig.[^ We apply 
post-selection and accept the output of the PBS only if 
there is a photon in both spatial output modes. 

Consider the effect of the PBS acting on a product 
state of two incoming photons: 

(a \H) + /3 |P))i 0 (a \H) + P\V ))2 ^ 

\H)^ \H)^ + |l/)3 |P)4 + aP{\H, V)^ + \H, V),) 

After post-selection the remaining quantum state is 

N{a^\H),\H)^+P^\V),\V),) = 

^ ((a^ \H)+P^ |P))3 |+)4 + (a" \H)-P^ |P))3 l-)^ 


where |±) = {\H) ± \V))/'/2 and N = 1/ is 
a norming factor. The success probability of the protocol 
is 1/A^ > 1/2. If we measure the photon at output mode 
4 the other photon collapses to N ■ (a^ \H) ± jP)) cor¬ 
responding to the measurement result. To get a def¬ 
inite outcome \H) + P"^ \ V) we may apply a Pauli- 
Z gate whenever we measure |—) (Fig. |l(f)[ ) or sim¬ 
ply neglect suc h case s introducing another level of post¬ 
selection (Fig. |l(g)[ ). Either way the protocol imple¬ 
ments a non-linear transformation S : a \H) P\V) i—>■ 
Npa^ \H)+P'^ |P)) which maps the identical qubit states 
of an ensemble to another qubit state of a smaller identi¬ 
cal ensemble. If we iterate this process S amended with 
an additional unitary step [/, the iterates (US)^ exhibit 
increasingly rich dynamics. 

It was shown m that iteratively applying US on an 
ensemble of identically prepared qubits the one qubit 
state of the ensemble \tpn) = ctn \H) + Pn \V) after n 
iterations may evolve sensitively with respect to the ini¬ 
tial state IV'o)- Nonetheless this was only shown for 
initial states lying on a zero-measure fractal on the 
Bloch sphere called the Julia set, see Fig. How¬ 
ever certain choices of unitaries from the family = 
cos(0) sin(d)e®‘^ 


— sm(0)e cos(0) 


seem to produce increasingly 


saturated Julia sets suggesting that it may reach a point 
where the whole Bloch sphere is covered by unstable ini¬ 
tial states. A candidate for such a transformation is 

[l i' 

^ = ^^7rl4,7r/2S where C/^/4,,r/2 = ^ 


Exponentional sensitivity for all initial states. In 

order to handle the arising non-linear maps better we 
project the surface of the Bloch sphere to the complex 
plane using stereographic projection. Thus a (photonic) 
qubit |'0) = a \H)+P \V) may be described using a single 
complex parameter z = ^ gC = CU{oo} including 
infinity. This representation yields (~^) a new description 
of our protocol in terms of rational functions mus]: 


S s : z —>■ z'^; ^ f : z 


z'^ + i 
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(d) e = 1 . 05 , 


f (e) 61 = 0.5, ¥> = 0.5 (f) 61= 0.232, ¥>= 0 


Figure 2: Julia sets consisting of the unstable initial 
states lying on the Bloch sphere. Julia sets (blue) are 
plotted for various choices of Ug^^p. (a)||(c) Shows how the Ju¬ 
lia set starts co verin g the whole Bloch sphere while increasing 
the value of 9. (d) Transition from a web to a simple closed 
curve, (e) The unstable states form a completely disconnected 
set. m The result of a parabolic explosion (implosion) |2U| 
where a stable fixpoint has become unstable - breaking a circle 
like connected Julia set into infinitely many parts. 


Using this formalism it turns out that / is one of 
a few special so called Lattes maps m and as such 
gained a lot of attention in the theory of complex dy¬ 
namical systems m- We can better understand the spe¬ 
cial properties of our Lattes map by analysing its rela¬ 
tionship to the corresponding linear transformation of 
the 2 dimensional torus. We will represent the torus 
C/Z[i] as the complex plane modulo the Gaussian in¬ 
tegers Z[i] = {a + bi\a,b G Z}. Its transformation is 
represented by multiplication with •(! — i): z —>■ (1 — z) • z 
mod (l,z) which rotates and folds the torus 2 times over 
itself. The correspondence between the torus and the 
sphere is established via the so called Weierstrass elliptic 
function p : C/Z[z] i—)■ C [55]. Relating the two surfaces 
gives rise to the identity = p((l — z)” • p~^{z)) 

showing that iterating / on C has essentially the same 
effect as repeatedly applying multiplication •(! — z) on 
C/Z[f], see Fig.[| 

Viewing $ through these glasses it becomes clear that 
it shows chaotic behaviour on the whole Bloch sphere. 
The map representing $ on the torus uniformly stretches 
the surface of the torus by a factor of and folds over 
itself two times. It is intuitively clear that the itera¬ 
tive application of such a transformation shows exponen¬ 
tial sensitivity to the initial position on the torus and 
has a positive Lyapunov exponent. Even more strikingly 
it exhibits exponential mixing, yielding that close ini¬ 
tial states separate exponentially fast on the surface of 
the Bloch sphere, as depicted on Fig. [^ For a rigorous 
derivation of the exponential mixing see Appendix. 



(e) ir^z)i > 1 (f) \rHz)\> 1 (g) ir«(z)i > i (h) \r\z)\> i 



(i) ir*(^)i>i (j) ir®wi>i (k)ir^°(z)i>i(i) ir“(z)i>i 


C/Z[z] : ^ ^ • • • -^ 




Bloch sphere : O -^ O 


p 

p 

p 

^ V P 






(m) A commutative diagram explaining the action of 3>. 


Figure 3 : Iterations of an exponentially mixing map. 

(a) ^ Visualisation of the iteratives of /. Each subfigure 
shows z € [—2, 2] X [—2z, 2z]; the domains are coloured accord¬ 
ing to whether |/°"| > 1 (black) or < 1 (white) distinguishing 
the northern and southern half of the Bloch sphere. After a 
few iterations even very close states get mapped to different 
halves of the Bloch sphere as indicated by the rapid alter¬ 
nation of black and white domains, (m) The iterative maps 
>1>°" act on the Bloch sphere correspondingly to the action of 
multiplication by (1 — z)" on the torus, explaining the regular 
pattern. 


Emergence of general complex rational dyanmics 
and the Mandelbrot set. We may find even more 
exotic transformations by generalising the protocol al¬ 
lowing us the physical realisation of any rational map 

z —?> of degree n > 2. The generalised scheme 

proceed by forming n-tuples 0 ”^i(q! |0) -I-/3 |l))j of iden¬ 
tical pure qubits and applying an appropriate n qubit 
unitary V. The final step is measuring all the qubits ex¬ 
cept the last one and keeping it only if all measurements 
resulted 0; this post-selective step can be shortly de¬ 
scribed by the projection [|0)(0 ... 00| -h |1)(0 ... 01|]. Im¬ 
plementing a specific rational map reduces to finding a 
V unitary satisfying: 

[|0)(0... 00| + |1)(0... 01|] U (G”^i(a |0) + /311)),) 

C n n \ 

Y, Cfco"/?"-" |0) + ^ |1) 

k=0 k=0 / 
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with arbitrary 7 7 ^ 0 , providing the desired post selected 
state. The existence of such V follows from a simple 
linear algebraic argument explicated in Appendix. 

A notable consequence is that we found a direct quan¬ 
tum physical realisation of the Mandelbrot maps 2 ; —>■ 
z'^ + c and can devise a quantum circuit for it. Fig. 
shows a possible quantum circuit implementation for this 
family of maps. The scheme demonstrates how the corre¬ 
sponding family of 2 qubit unitaries may be constructed 
using only controlled 1 qubit gates which are considered 
experimentally more feasible in general. 


1 ^) 

1 ^) 


H 




-& 


(n]— -ip —\x 


-^=r = 0 


Figure 4: A quantum circuit implementing Man¬ 
delbrot maps. The controlled rotation of the gate la¬ 


belled by Tj is 
belled by (p is 


I 0 ■ 
0 


1 


-1 


and the phase gate la- 


while H stands for a Hadamard 


gate. If we set ri = l/r 2 , r 2 = |c| • ^{I + ^ + 4/|cp)/2 
and then the resulting map is z^ + c^ if we accept 

only the 0 measurement outcome on the second qubit. 


The cost of non-linearity and the Schrodinger mi¬ 
croscope. We have just shown that using post-selection 
one can implement a wide range of non-linear maps that 
may be useful for various tasks. A highly non-linear 
map like we presented provides a sort of ’’Schrodinger 
microscope” [53] enabling us to exponentially magnify 
tiny differences between quantum states. Having such a 
tool we are tempted to develop powerful quantum al¬ 
gorithms utilising it. It is well known that introduc¬ 
ing post-selection to Quantum Computing makes it ex¬ 
tremely powerful - the corresponding complexity class 
PostBQP [24] includes NP and even PP. The question 
of efficiency and resource needs naturally arises. We ad¬ 
dress it using a black box argument considering results 
about state discrimination [531157] . 

Suppose we have a quantum device implementing n 
iterations of d* processing a qubit ensemble of size N. 
The size M of the successfully processed output ensem¬ 
ble may be probabilistic. We would like to determine 
its success rate i.e. the average ratio M/N. To de¬ 
rive a bound on the success rate consider applying this 
quantum device to the qubit ensemble having state ei¬ 
ther 1 ^ 0 ) or I'i/ji) 7 ^ 1 ^ 0 ) with equal probability 1 / 2 . 
If their distance d(|V'o) , l^i)) = \/l - | (V’olV’i) P < 
1 /A^ then the distance of the full ensemble states is 

d({g)^ IV’o), 0^ IV’i)) ~ VN ■ d(IV'o), ll/i)) thus we can¬ 
not distinguish the possible ensemble states with error 
probability less than « (1 — VN ■ dd'f/o), |V’i)))/2, see 
[551 Chapter IV §2]. Suppose d{\tjjj) ,\+)) 1/2” for 

j = 0,1 then, after n iterations of our process the dis¬ 


tance of the two sates increases by a factor of roughly 2 ” 
since $ doubles infinitesimal distances around the fixed 
point |-|-). If the device outputs M copies of the trans¬ 
formed states then we can distinguish the ensembles with 
error probability « (1 — '/M ■ 2” • ddt/o) j l'0i)))/2. For 
large N the success rate is roughly constant, because of 
the law of large numbers, thus we can treat the value 
M/N fixed. But y/M ■ 2"/-\/]V cannot exceed 1 as the 
error probability of discrimination cannot decrease and 
so the success rate is upper bounded by 4“”. This holds 
for states lying close to |-|-), in better cases the rate may 
be higher. For our implementation scheme each iteration 
has a success rate at least 1/4 (up to a negligible term 

— 1 /N due to parity) implying that this scheme provides 
the best possible worst case success rate. 

In this way we showed that exponentially many copies 
are needed for n iterations of the process. Similar upper 
bound can be devised to any non-linear map that have 
a region where the separation of close states can be de¬ 
scribed by a multiplicative factor A > 1. If we follow the 
above argument it turns out that the worst case success 
rate of such protocol is bounded by 1/A^. Note that the 
particular choice of metric by which we measured separa¬ 
tion is not limiting the scope of the argument too much 

- we could use any other metric which agrees infinitesi¬ 
mally, e.g. the Bures metric. Thus it turns out that the 
implementation of any kind of Schrodinger microscope 
needs exponentially many copies of the states in terms of 
magnification steps, more precisely quadratically many 
in terms of the total magnification. 

Discussion 

While exploring the possible dynamical properties of 
state selective protocols we found that any complex ratio¬ 
nal map can be implemented using state selection. Such 
a general and natural correspondence between a physi¬ 
cal system and the theory of complex dynamical systems 
is unique up to our knowledge. We could also devise a 
realistic optical experimental scheme which implements 
particularly interesting quadratic rational dynamics. 

Meanwhile we studied the emergence of exponential 
sensitivity and it lead us to the analysis of implementa¬ 
tion cost which turned out to be exponential. We found 
a general bound on the number of copies needed for the 
successful implementation of any expanding non-linear 
map. We proved that a protocol capable of magnify¬ 
ing differences between close quantum states by a factor 
A > 1 necessarily yields a rate of loss at least 1 — 1/A^ 
in the number of copies of the unknown input quantum 
states. This ” Quantum magnification bound” is basically 
another reformulation of the fact that one cannot boot¬ 
strap quantum information without an external source 
and somewhat resembles the quantum no-cloning theo¬ 
rem. 

We used the ’’Quantum magnification bound” prin¬ 
ciple to show the optimality of our implementation of 
an exponentially mixing map. Similarly, it may be ap¬ 
plied to other quantum information protocols providing 
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a general tool for bounding the success rate of particu¬ 
lar probabilistic protocols. This principle also helps to 
understand the advances and limitations of any kind of 
Schrodinger microscope regardless the actual implemen¬ 
tation method. The term Schrodinger microscope was 
first introduced for a protocol using collective weak mea¬ 
surements and coherent feedback |23] where the limita¬ 
tions were unclear due to the approximate arguments 
applied. A fundamental approach like the one we pre¬ 
sented should be helpful for understanding such complex 
systems as well. 

Looking at general processes with inspiration coming 
from this principle may also provide a new insight to the 
relation of classical and quantum chaos [28], suggesting 
that classical deterministic chaos may be just an approx¬ 
imation with a characteristic time scale. Classical deter¬ 
ministic chaotic systems explode fast and observing the 
deterministic evolution of the system even at a macro¬ 


scopic level enables the determination of the initial con¬ 
ditions increasingly precisely [29] . But there is a level of 
precision that is prohibited by quantum uncertainty re¬ 
lations. This is an apparent philosophical contradiction 
if we believe classical physics is based on quantum me¬ 
chanics; a possible dissolution is saying that on long time 
scales one cannot treat a classical process deterministi¬ 
cally chaotic just chaotic in some statistical sense. 
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Appendix 

2 . - 

Connection to Lattes maps. The map / = can 

be written as o o , where / = 3 ^ is a well- 
known Lattes map [ 2 T] , and is a self inverse Mobius 
transformation [19]. Since conjugation by Mobius trans¬ 
formation do not change the iterative features it implies 
that / exhibits the same dynamics as /. To give a phys¬ 
ical meaning to this Mobius transformation we mention 
that it corresponds to a rotation of the Bloch sphere i.e. 
/ and / essentially describe the same process just written 
in a different qubits basis. 

As we already indicated in this article / is conjugate 
to the map (1 — i) • idc/z[i] via the Weierstrass-p func¬ 
tion. In fact, we need a slightly transformed version of 
the Weierstrass-p which amended by our Mobius trans¬ 
formation. Our transformed version can be written as 
p(z) := o definitions of L,g 2 and 

p(z; L) see jH]). Using this function we get the key iden¬ 
tity /°” O p = p o (1 - i)" • idc/z[i] ■ 

p induces a two sheet branched covering with 4 ex¬ 
ceptional points which are covered only once. (These 
points can be easily spotted on Fig. |^. With the excep¬ 
tion of these 4 points p is a two-to-one map so it does 
not have a well defined inverse. A general point i S C 
has pre-images z, —z € C/Z[f], but the linear map on 
C/Z[z] carries opposite numbers to opposite ones so for¬ 
tunately the identity /(z) = p((l — i) ■ p“^(z)) Vz G C 
holds regardless which branch of p~^ is considered. In 
this sense we can say that the stronger identity f°^{z) = 
p((l — f)” • p“^(z)) Vz G C also holds. 

Metric on the Bloch sphere C and the torus C/Z[f]. 
We would like to show exponential mixing of our Lattes 
map thus we need to understand how distances are dis¬ 
torted by p. In order to trace the problem we need to 
introduce some proper distance concepts. 

A possible metric on pure quantum states is given 
by using the distance defined by the quantum angle 
dA{\'4’zo) I \4’z)) '■= arccos | {ipoli’z) \- Note that this dis¬ 
tance coincides with the natural spherical metric of the 
Bloch sphere up to a multiplicative factor of 2. This 
metric is similar to the Bures metric defined by the 
distance ds = , |V' 2 )(V'z|)), where 

F{.,.) is the Fidelity of two density matrices. A third 
possible distance definition is given by d{\tpzo) > IV'i)) = 
\/l — I (^zolV'z) P- Since all three distance definitions co¬ 


incide for infinitesimal distances we are free to chose the 
most appropriate one for our calculations. For now we 
stick with the natural metric of C which is 2 times the 
quantum angle ^^(zo, z) = 2 • dA{\'4’zo) j the index 

R refers to the fact that dn is a Riemannian metric. For 
the torus C/Z[i] we use the natural Riemannian metric 
inherited from C. 

First we show that the Mobius transformation 
leaves the metric on C invariant. I.e. we would like 
to show that the conformal metric p := is triv¬ 

ial, where ds and ds' are tangent vectors of C at 
points s and s' = such that ds is mapped to ds' 
by the tangent map. We proceed using the identity 
~ ^here dz and dz' are tangent vec¬ 

tors of C. The conformal metric transformation in¬ 
troduced by the stereographic projection is well known 
^ similarly . Finally 


to be - _ 

\dz'\ — |z'|2-|-l’ 

d[^)ldz 


|ds| 

(z-l)-(z-H) 
(z-l) = 


2 

|z-l| = 


Putting 


|ds| 4 |z-l|2 |fL[|2 + l 

= 1 as we indicated. 


everything together J|^ = 

hl^+i 2 |z-i|^ 

1 |z-l|2 p+l|2 + |^_l|2 

The Mobius transformation is just an isometry of C 

so we can concentrate on the other part of 

*V92(i'=*) 

our Weierstrass function p. Now p := 

where dz' is a tangent vector of C at a point z'. Just as 

above ^ = ^ 7 ^ = ^|p(t)| 2 +i - For the other factor 


W = 

well known property p'^ = 4p^ — < 72 P — < 73 - Since (73 (L = 
t) = 0 the final formula is: p^ = ^ 

g 2 p{z) \ = 64g2 '(I|p(l) 72 +g 2 p' ■ Using the triangle inequal- 

ity we get p < 64g2 ( 4 |p(z)| 2 +g 2 )= “ 64(72■ 

This function has its maximum when |p(z)p = 52 / 4 , 
substitution yields p^ < 16yd72. Finally using that 
y/g 2 ~ 13.7504 we arrive at the conclusion p < 16. 

The conformal metric p < 16 is upper bounded mean¬ 
ing that the image of any two points from the surface of 
the torus gets mapped to points having spherical distance 
less than 16 times their torical distance. Looking at this 
inequality the other way round shows that if have a point 
s on the surface of the sphere and a radius s ball (s, e) 
around it then this ball’s pre-image p“^(iJdj,(s, e)) con¬ 
tains another ball Bd^^^^.^{p~^[s),e/1&). 


Exponential mixing. Suppose we have a one qubit 
state \ipz) prepared with accuracy e meaning that it is 
a pure quantum state \'tpz) close to some jV’zo) such that 
the quantum angle dAdV’zo) > IV’z)) = arccos | (V’zolV’z) I < 
£. In other words |'0z) lies in the diameter e ball 
around \ilizo) i-e- IV’z) € (IV’zo) > This ball cor¬ 
responds to Bdj^{zo,2e) using the (Bloch) spherical rep¬ 
resentation. As we already discussed p“^(i3d^(s, 2 e)) D 
Udc/z[.](p”^(s),e/8). 
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It is easy to show that after n = — log^(e/8) iteration 
of •(! — i) over all the points of (p“^(zo), <5) covers 

the whole C/Z[i]. Thus for n = 6 — log^(e) we have 
C = f°^{Bdj^{zo,e)). It means that after 6 — log^(e) 
iterations the initial uncertainty about the state |'0z) 
evolves so much that $°"(|^/>z)) may be any pure state. 
This statement is basically a translation of the fact that 
the linear map -(I —i) has Lyapunov exponent ln(-\/2) on 
C/Z[z] and shows exponential mixing and sensitivity on 
the surface of the Bloch sphere. 

Construction of n-qubit unitaries for degree n ra¬ 
tional maps. We would like to implement the rational 

CkZ^ 

function z —>■ ■ Our generalised protocol starts 

by forming n-tuples of identical pure qubits of our en¬ 
semble then continues by the application of a specific n 
qubit unitary V. The hnal step is a measurement on all 
the qubits except the last one of every tuple. The pro¬ 
tocol succeeds if all the measurements resulted in 0, the 
unmeasured qubit is kept only in such cases. 

Initially the state of the n-tuples is the following prod¬ 
uct state: 


l<^) = (g)(« |0), + /? |1),) = E \4>k) 

i—1 k—0 

where \(t)k) = \b) 

bG{0,l}" 

As before we use the parametrisation z = ajfi for a 
qubit (o |0) -I- /3 |1)). Then the parameter of the unmea¬ 
sured, post selected qubit can be described as follows: 


n—1 zeros Cfc 

( r'.'o 0|C|<^) 

1 ■v' (1) 

n—1 zeros ^ 

ELo4a'=/3"-'= ELo^^'^ 

We present a linear algebraic argument showing that 
for any rational map of degree > 2 there is a suitable uni¬ 
tary V, thus we need to find a unitary for any coefficients 
Ck,dk describing a rational map. 


First let us introduce some vector jw) = (| 10) — 

n—2 zeros 

I 01))/-\/2 which is orthogonal to all the \(j)k) vec- 

n—2 zeros 

tors. We set vq = J2k=o -I- a; • |a;) and vi = 




a/O) 

+ y ■ |w). Then by choosing x,y G C ap- 


propriately we can always satisfy the equalities: 


2/*a; = - ^ cldk ^ cldk +y*x = 0-^ (do|Fi) = 0 
k^O fc=0 

( 2 ) 

n 

n n ^ 

^ E = E 

k^O 

Finally setting ( 0| ^ •= (^o| / ('^^oKo) 

n—1 zeros 

( 1| y ■= (All / (Ai|wi) satishes Q and also ©,(§ 

n—1 zeros 

so we can extend F to a full n qubit unitary by dehning 
the remaining 2" — 2 orthonormal rows arbitrarily. 

Note that if EEo EEo common 

roots, than the probability that the process succeeds is 
greater than some p probability regardless the state z. 
The probability cannot be zero since there is no common 
root, and it follows from compactness that is cannot ap¬ 
proach zero for a fixed map. However, depending on the 
map, this lower bound may be arbitrarily low. 

Thus using the above defined V unitary we can imple- 

E”’ CkZ^ 

ment the rational function ^ where as before this 

means a transformation 


(o|0)+/3|l)) 


N E |0) -f E 40'=/?"-'= |1) 






where z = ajp and is a norming factor. 














